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Abstract 

The  Generalized  Implicit  Flow  Solver  (GIFS) 
computer  program  has  been  modified  and  applied 
for  analysis  of  three-dimensional,  reacting,  two- 
phase  flow  simulation  problems.  The  intent  of  the 
original  GIFS  development  effort  was  to  provide 
the  Joint  Army,  Navy,  NASA,  Air  Force  (JANNAF) 
community  with  a  standard  computational  method¬ 
ology  to  simulate  multiple  nozzie/plume  flow-field 
phenomena  and  other  three-dimensional  effects. 
The  Van  Leer  Flux  Splitting  option  has  been  suc¬ 
cessfully  implemented  into  the  existing  GIFS 
model  and  provides  a  more  robust  solution 
scheme,  making  application  of  the  model  more 
reasonable  for  engineering  applications. 

This  paper  is  a  continuation  of  the  previous  work 
and  reports  the  significant  results  of  parametric 
flow-field  simulations  resulting  from  a  dual-nozzle 
propulsion  system  operating  at  a  high-altitude  flight 
condition.  In  support  of  this  effort,  four  calculations 
of  Titan  II  SLV  flow  fields  have  been  completed  to 
assess  the  effects  of  three-dimensionality,  missile 
body  effects,  chemistry,  and  gas  generator  flow  on 
simulated  plume  exhaust  flow-field  properties. 

These  calculations  indicate  that  three-dimen¬ 
sionality  is  always  an  important  factor  and  could 
substantially  influence  the  interpretation  of  the 
results.  If  three  dimensional  effects  are  oversimpli¬ 
fied  in  the  model,  analyses  of  the  spatial  results 
can  be  misinterpreted  and  misapplied.  The  missile 
body  effect  can  also  generate  three-dimensional 
influences  affecting  the  Mach  reflection  location, 
the  plume/plume  impingement  shock  location, 
inviscid  shock  structure,  and  shear  layer  growth. 
The  gas  generator  flow  influences  the  very  near¬ 
field  simulation  but  is  significantly  dampened 


beyond  approximately  1  meter  downstream  of  the 
nozzle  exit  (nonreacting).  If  missile  base  heating, 
recirculation  effects,  or  nozzle  impingement  heat- 
transfer  analyses  are  required,  the  gas  generator 
flow,  including  a  kinetic  chemistry  model  is  the 
dominant  influence. 

introduction 

In  order  to  support  propulsion  testing  and  analy¬ 
sis  requirements  of  the  aerospace  and  exhaust 
plume  community,  a  need  exists  for  a  fluid  dynam¬ 
ics  model  that  solves  the  fully  coupled  two-phase 
Navier-Stokes  equations  in  multiple  dimensions. 
Evaluations  of  solid-propellant  rocket  motor  perfor¬ 
mance,  nozzle  erosion,  and  rocket  plume  radiative 
transfer  analyses  require  a  computer  model  that 
simulates  complex  three-dimensional,  chemically 
reacting,  two-phase  flow  effects.^  Although  this 
type  of  full  Navier-Stokes  method  provides  an 
accurate,  qualitative  description  of  the  basic  fea¬ 
tures  of  the  propulsion-generated  flow  fields,  quan¬ 
titative  simulations  for  predicting  fundamental 
parameters  such  as  base  pressure  and  heat  trans¬ 
fer,  gas  static  pressure,  temperature,  and  chemical 
composition  in  the  flow-field  domain  have  not  been 
validated.  In  the  past  few  years,  significant 
progress  has  been  made  in  the  areas  of  numerical 
rocket  flow  simulations  and  computational 
resources  to  the  point  that  Navier-Stokes  solutions 
are  viable  analysis  tools.  Efforts  to  validate  these 
models  tend  to  lag  the  algorithm  development  and 
the  computational  hardware  advancements  consid¬ 
erably  and  hamper  the  overall  confidence  of  the 
quantitative  results.  This  study  provides  a  qualita¬ 
tive  validation  of  the  physical  effects  included  in  the 
GIFS  model  and  provides  a  basis  for  in-depth 
understanding  of  the  results  of  the  GIFS  model. 
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The  flow  fields  generated  by  rocket  propulsion 
systems  are  compiex,  with  regions  of  strong  invis- 
cid/viscous  interactions,  free-stream  shear  layers, 
nozzle  wall  and  missile  body  boundary  layers, 
external  and  internal  shocks,  separation  regions, 
and  plume/plume  impingement  and  associated  flow 
interactions  for  multiple  nozzle  designs,  all  with 
chemically  reacting  kinetics.^  Coupling  all  of  these 
phenomena  simultaneously  in  a  numerical  simula¬ 
tion  tool  challenges  the  state  of  the  art  (SOA)  for 
CFD  models.  To  account  for  all  phenomena  affect¬ 
ing  plume  flow  properties  and  the  resulting  radiative 
transfer  implications,  computationally  efficient,  mul¬ 
tidimensional  computer  models  are  required. 

Recenty  there  has  been  increased  emphasis  on 
the  application  of  CFD  models,  both  in  the  com¬ 
mercial  arena  and  government-developed  com¬ 
puter  programs,  to  simulate  complex  multidimen¬ 
sional  plume  phenomena.  A  large  majority  of  the 
solutions  obtained  to  date  are  based  on  the  perfect 
gas  (constant  gamma)  approximation,  as  fully 
reacting  flows  are  considerably  more  complex  and 
difficult  to  solve.  Including  the  effects  of  chemistry 
in  the  solution  produces  a  stiff  set  of  equations  that 
are  numerically  difficult  to  solve  using  conventional 
algorithms.  In  addition,  the  grid  resolution  require¬ 
ments  become  more  severe,  and  time  step  issues 
arise  when  reacting  chemistry  is  included  in  the 
Navier-Stokes  model. 

Conventional  rocket  exhaust  flow  computer 
models^’"*  employ  Euler  solutions  for  the  plume 
core  flow  and  superimpose  a  turbulent  parabolic 
mixing  approach  for  the  free-stream  air  and  plume 
entrainment  (shear)  layer.  These  methods  are  gen¬ 
erally  not  adequate  for  situations  when  the  flow  is 
not  fully  dominated  by  either  inviscid  core  expan¬ 
sion  or  the  plume  afterburning  phenomena.  Three- 
dimensional  features  produced  by  multiple  nozzle 
propulsion  systems  and  vehicle  body/base  interac¬ 
tions  cannot  be  sufficiently  treated  through  the  use 
of  these  models.  Inaccurate  accounting  of  the  3D 
upstream  influences  on  the  plume  shear  layer 
development  and  the  resulting  plume  structure  is 
one  of  the  primary  issues  preventing  approximate 
models  from  accurately  simulating  the  overall  flow- 
field  phenomena.^ 

The  GIFS  numerical  algorithm  provides  a  solu¬ 
tion  of  the  two-  and  three-dimensional  Reynolds- 


averaged  Navier-Stokes  (NS)  equations  using  the 
MacCormack  implicit  finite-volume  algorithm  with 
Gauss-Seidel  line  relaxation.®  Several  2-D  and  3-D 
plume  flow-field  calculations  have  been  completed 
for  the  plume  near-field  region  using  the  original 
version  of  the  GIFS  model.®  The  GIFS  model 
includes  a  frozen  and  generalized  finite-rate  kinetic 
chemistry  model,  a  Lagrangian  particle  model  for 
treating  solid  and  liquid  particulates,  and  a  two- 
equation  turbulence  model,  as  well  as  a  laminar 
model.  These  complex  phenomena  are  required  to 
accurately  simulate  the  physics  expected  to  con¬ 
tribute  to  the  flow-field  spatial  distribution  of  gas 
dynamic,  thermodynamic,  and  chemicai  properties. 
The  Van  Leer  Flux  Splitting  option^  has  been  suc¬ 
cessfully  implemented  into  the  original  GIFS  model 
and  provides  a  more  robust  solution  scheme  for 
simulating  propulsion  flow-field  phenomena.®  This 
enhanced  version  of  the  GIFS  model  was  applied 
in  this  study. 

The  near-field  plume  flow  field  emanating  from 
a  multinozzle  vehicle  flying  at  high  altitude  are 
largely  dominated  by  the  processes  taking  place  in 
the  plumes’  interaction  region.  Prior  to  the  release 
of  the  GIFS  model,  simulations  of  multiple  nozzle/ 
plume  flow  fields  were  commonly  treated  by 
assuming  a  single  equivalent  nozzle  configuration 
having  equal  mass,  energy,  and  momentum  of  the 
multiple  nozzle  geometry.  Further,  uniform  (one¬ 
dimensional)  nozzle  exit  flow  properties  were  used 
as  the  starting  conditions  for  the  plume  calculation 
as  a  simplifying  assumption.  The  simplified  model 
assumes  that  the  details  of  the  3-D  flow  structure  in 
the  near-field  flow  and  the  2-D  start  condition  are 
unimportant  and  that  the  flow  processes  affecting 
the  plume  shear  layer  initialization  (such  as  base 
separation  and  recirculation)  will  be  dominated  by 
the  overall  ambient  flow  entrainment  effects.  For 
spatial  analysis  requirements  the  level  of  agree¬ 
ment  between  computations  based  on  the  single 
equivalent  nozzle  methodology  and  simulations 
from  multiple  nozzle  propulsion  systems  has  not 
been  acceptable.  The  source  of  the  disagreement 
is  due,  at  least  in  part,  to  an  incorrect  physical 
model  of  the  phenomena  dominating  the  observa¬ 
tions,  e.g.,  inadequate  turbulence,  incomplete 
chemical  mechanisms,  missing  or  inaccurate  reac¬ 
tion  rates,  simplified  initial  start  conditions  and 
three-dimensional  geometry  effects.  The  actual 
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physical  geometries  are  often  oversimplified,  and 
relevant  details  of  the  enQine  design  and  the  result" 
ing  impact  on  the  initial  conditions  and  the  down* 

Oilmen  1 1  OAiictuoi  iiv/vv  ai^  iicv^uoiuiy  i^iiwiou. 


In  the  previous  study, ^  six  simulations  of  the 
Titan  II  Plume  flow  field  were^  completed  to  assess 
the  effects  of  three  dimensionality, 


cosity,  finite-rate  chemical  kinetics,  intranozzle 
geometric  spacing,  and  nozzle  exit  profile  initial 
conditions.  This  study  is  a  continuation  of  the  previ¬ 
ous  work  and  further  demonstrates  the  signiiicance 
of  three-dimensional  effects  as  applied  to  multiple- 

T"u.«. 

luurvot  iiiiooiio  piuiiie  oiiiiuiauui lo.  i  iio 

results  can  be  applied  to  guide  efforts  to  determine 
how  these  phenomena  may  be  approximated  to 
improve  engineering  approaches,  and  to  explore 
how  computer  resource  regulrements  iwi  oniiuioi- 
mg  three-dimensional  solutions  can  be  reduced. 
The  results  presented  in  „  „  „  -  nu 
this  paper  demonstrate 
continued  progress  in 
simulating  the  rocket 
exhaust  flow  field  from  a 
multiple-nozzle  propul¬ 
sion  vehicle  in  flight. 


zen  flow  conditions.  To  assess  the  effects  of  the 
g35  Q8n6rat0r,  tho  SoCOnci  C358  WaS  r8p8at8Cl  Wiih" 

out  consideration  of  the  gas  generator  in  the 

liiv/uct.  I  no  111101  owiuiioii  ooouMiovi  o  Oil  1^10  OlyfUIVO- 

lent  axisymmetric  nozzle  to  approximate  the  three- 
diiTiensional  solution  geometry.  Included  in  Table  1 
is  a  summary  of  all  the  parametric  flow-field  condi¬ 
tions  and  assun'iptioi'is  evaluated  in  this  study  and 
the  previous  investigation.  Cases  1-6  were 
reported  in  the  previous  paper."  This  study  focused 
on  the  last  four  cases.^"’® 


The  Titan  11  Space-Launched  Vehicle  (SLV) 

|jiu|^uioiuii  oyoioiii  to  ^lujjciiou  uy  ivvu  on^iiioo 

located  on  either  side  of  a  plane  of  symmetry  pass- 
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lil^  UIIOU^II  II  10  VOIII\^IO  UOI  KOI  III  10.  \sillVOII  III  1  1^.  I  10 

a  schematic  of  the  vehicle.  The  sea-level  rated 
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Computations 

This  computational 
effort  consisted  of  four, 
three-dimensional  twin- 
nozzle  calculations  for 
the  Titan  II  vehicle  at 
actual  flight  conditions. 

The  2D-axisymmetric 
initial  conditions  were 
computed  using  the  Two-Dimensional  Com¬ 
puter  Program,  TDK  (Ref.  9)  and  were  identi¬ 
cal  for  all  cases.  For  all  cases,  the  plume  flow 
field  was  simulated  at  an  altitude  of  approxi¬ 
mately  50  km  at  a  free-stream  Mach  number 
of  5.7,  assuming  turbulent  conditions  in  all 
flow  regions.  The  first  case  simulates  three- 
dimensional,  reacting  flow  conditions  and 
includes  flow  regions  surrounding  the  missile 
body,  the  base  flow  region,  the  gas  generator 
flow,  and  the  exhaust  plume  domain  assum¬ 
ing  reacting  chemistry  conditions.  The  sec¬ 
ond  case  is  similar  to  the  first,  assuming  fro¬ 


0.8  Rb 


Nozzle  Cant  Angle  =  2  deg,  1  min 


Engine  Gimbal  Point 

Fig.  1.  Vehicle  dimensions. 

Table  1 .  Parametric  Flow-field  Conditions 


Case 

Viscosity 

iniet 

Profile 

Chemistry 

Nozzle  Spacing 

GG 

1 

HgjjQjgg 

1-D 

no 

2 

1-D 

Constant  y 

wide 

no 

3 

Turbulent 

1-D 

Constant  j 

Narrow 

no 

4 

Turbulent 

1-D 

Finite  Rate 

Narrow 

no 

5 

Turbulent 

1-D 

Narrow 

no 

6 

Turbulent 

2-D 

Narrow 

no 

mm 

Turbulent 

Finite  Rate 

Narrow 

8 

Turbulent 

Frozen 

Narrow 

wm 

9 

Turbulent 

2-D 

Frozen 

Narrow 

no 

10 

Turbulent 

2-D 

Finite  Rate 

2D-axisymmetric 

no 

vided  by  gimballing  of  the  engine  5  deg  from  the 
engine  neutrs!  position  (2“deg  cent  angle  away 
from  the  vehicle  centerline  for  both  nozzles).  The 


plume  configuration  is  shown  in  Fig.  2b.  The  exter¬ 
nal  airflow  conditions  and  nominal  1-D  liguid-pro- 
pellant  rocket  nozzle  throat  and  exit  conditions  and 


2-deg  cant  angle  was 
included  in  the  three- 
dimensional  simulations.^® 
Both  nozzles  are  identical 
in  geometry  and  operating 
conditions. 

The  thrust  chamber 
assemblies  and  the  nozzle 
skirt  are  regeneratively 
cooled.  In  addition,  injector 
spray  patterns  intention¬ 
ally  direct  a  fuel-rich  layer 
adjacent  to  the  chamber 
walls  to  reduce  the  wall 
heat  loads.  Power  to  drive 
the  turbopumps  is  derived 
from  two  gas  generators 
which  represent  approxi¬ 
mately  1.5  percent  of  the 
overall  propellant 

expended  by  both  engines. 
These  gas  generators  are 


intentionally  operated  at  a  Schematic  of  the  3-D  grid  and  the  planes  of  symmetry 

fuel-rich  conditions  to  mini¬ 


mize  heat  loads  on  the  tur¬ 
bine  blades. 

To  limit  the  number  of 
grid  points  and  minimize  the 
CPU  execution  time,  quar¬ 
ter-plane  symmetry  assump¬ 
tions  were  made  for  the  com¬ 
putational  domain.  All  simu¬ 
lations  were  accompiished  at 
zero  angle  of  attack.  A  total 
of  4  million  grid  points  was 
utilized  in  the  computational 
domain  for  the  three-dimen¬ 
sional  simulations,  and 
200,000  grid  points  for  the 
axisymmetric  case.  The 
exhaust  plume  portion  of  the 
computational  domain  is 
shown  in  Fig.  2a,  with  the  ori¬ 
entation  indicated  by  the 
axes.  A  schematic  of  the  3-D 


b.  Schematic  of  the  3-D  plume 
Fig.  2.  Geometry  representation  of  the  computation  domain. 
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gas  generator  conditions  for  the  Titan  cases  are 
presented  in  Table  2.  The  calculations  assuming 
finite-rate  chemical  kinetics  used  a  chemical  reac¬ 
tion  model  consisting  of  1 0  species  and  1 1  reac¬ 
tions  for  carbon,  hydrogen,  oxygen,  and  nitrogen- 
based  propellant  systems.  The  K-e  turbulence 
model  was  used  for  the  viscous  stress  approxima¬ 
tion.  The  computation  was  performed  on  a  single¬ 
processor  SGI  Power  Challenge  R8000.  For  the 
three-dimensional  cases,  the  calculations  required 
approximately  twelve  weeks  of  CPU  time  to  con¬ 
verge.  In  contrast,  the  axisymmetric  case  required 
5  days  to  converge.  The  GASP  computer 
program  11  contains  a  hybrid  Navier  Stokes/ 
Parabolized  Navier  Stokes  (PNS)  methodology 
which  was  also  applied  for  these  calculations.  The 
results  compared  favorably  with  the  GIFS  calcula¬ 
tions;  however,  the  computational  time  was  signifi¬ 
cantly  reduced  by  using  the  PNS  method  in  the 
supersonic  flow-field  regions.  For  these  cases,  the 
3D  calculation  time  was  reducted  to  four  weeks  and 
the  axisymmetric  to  18  hours  using  the  PNS  model. 

Table  2.  Inflow  Conditions 


Ambient  Conditions  at  47.6  km 

-  Tinf  =  269  K 

-  V=  1877.6  m/sec 

-  Cp/Cv=  1.4 

-  r=  1.388  X  10“^  kg/m^ 

-  P  =108  Pa 

-  Mach  =  5.7 

-  Species  Concentrations  (Mass  Fraction) 

N2  =  0.77 

O2  =  0.23 

Jet  Conditions  (One  Dimensional) 

-  T=1920K 

-  V  =  2776.6  m/sec 

-  P  =  92800  Pa 

-  r=  1.68  X  lO"”*  kg/m^ 

~  Mach  =  3.0 

-  Species  Concentrations  (Mass  Fraction) 

CO  =  0.039  GO2  =  0.1 81 1  H2O  =  0.3496 
N2  =  0.414  NO=  0.0109  OH  =  2.139e-3 

H2  =  3.13e-3  H  =  1 .24e-4  02=  0.0  O  =  0.0 


Gas  Generator  Conditions 


-  M  =  1.01 

-  T  =  899  K 

-  P  =  85488.9  Pa 

-  CO  =  0.04 

-  CO2  =  0.004 


-  CH4  =  0.135 

-  H2=  0.035 

-  H2O  =  0.034 

-  NH3  =  0.253 

-  C  =  0.038 


Including  the  previous  work,  ten  numerical  sim¬ 
ulations  of  the  Titan  11  twin-nozzle  configuration 


flow  field  were  obtained  at  Mach  5.7.  Different 
combinations  of  flow  assumptions  and  approxima¬ 
tions  were  applied  in  an  attempt  to  isolate  individ¬ 
ual  physical  influences  and  effects.  In  the  previous 
study,  only  the  exhaust  plume  flow  field  was  com¬ 
puted.  The  upstream  influence  of  the  missile  body, 
the  missile  base  region,  and  the  gas  generator  flow 
were  not  considered. 

The  numerous  assumptions  and  specifics  of  the 
all  the  solutions  are  summarized  briefly.  Cases  1 
and  2  incorporated  laminar  and  turbulent  viscous 
stress  models,  respectively,  and  were  computed 
assuming  a  perfect  gas  equation  of  state.  Case  3 
was  a  turbulent,  constant  gamma  approximation 
v^^ith  an  intranozzle  geometric  spacing  that  differed 
from  Cases  1  and  2.  Cases  4  and  5  differed  in  the 
chemistry  approximation,  assuming  frozen  and 
finite-rate  chemistry,  respectively.  Cases  1-5  all 
assumed  uniform  (1-D)  nozzle  exit  properties  as  the 
GIFS  start  line  conditions.  Case  6  assumed  a  2D 
start  line  profile.  Cases  1-  6  were  all  initiated  at  the 
nozzle  exit  location  and  did  not  include  upstream 
body/base  region  in  the  simulation.  The  nozzle  exit 
plane  (calculated  via  TDK)  was  used  to  define  the 
starting  boundary  condition  for  the  GIFS  plume  cal¬ 
culation.  Cases  7,  8,  and  9  include  upstream  influ¬ 
ences  and  differ  in  the  chemistry  assumption  and 
inclusion  of  the  gas  generator  in  the  simulation. 
Case  7  simulated  chemically  reacting  flow  over  the 
missile  body,  in  the  base  region,  and  the  exhaust 
plume  domains.  The  gas  generator  flow  was 
included  in  this  simulation.  Case  8  contrasts  the 
influence  of  the  chemistry  by  repeating  the  previous 
calculation  considering  a  frozen  flow  assumption. 
Case  9  is  similar  to  Case  8,  but  excludes  the  gas 
generator  flow.  Case  10  is  an  axisymmetric  single 
equivalent  nozzle  approximation  of  the  twin-nozzle 
geometry  and  includes  the  missile  body,  missile 
base  region,  and  exhaust  plume  in  the  computa¬ 
tional  domain.  The  gas  generator  flow  was  not  con¬ 
sidered  in  the  axisymmetric  simulation. 

The  following  will  highlight  the  effects  evident 
from  contrasting  the  various  solutions. 

Contrasting  the  turbulent  and  laminar  solutions 
indicated  that  the  barrel  shock  reflection  point  is 
located  approximately  8  meters  farther  down¬ 
stream  for  the  turbulent  solution.  Therefore,  the 
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spatial  characteristics  of  the  inviscid  plume  flow 
field  structure  can  be  influenced  by  the  viscous 
stress  model  included  in  the  solution.  In  the  present 
case,  at  47.6-km  altitude,  the  effect  of  turbulence  is 
not  particularly  strong.  However,  at  lower  altitudes 
where  higher  Reynolds  numbers  are  expected,  tur¬ 
bulent  effects  will  be  more  predominant. 

The  influence  of  intranozzle  spacing  was  also 
investigated.  The  results  indicate  that  intranozzle 
spacing  has  a  noticeable  influence  on  the  flow-field 
structure.  The  initial  plume  expansion  angle  is 
larger  for  the  wider  spacing  case,  and  the  shock 
reflection  location  is  farther  downstream.  It  is  con¬ 
cluded  that  the  intranozzle  spacing  affects  the 
plume  impingement  shock  location,  the  inviscid 
shock  structure,  and  the  shear  layer. 

The  effects  of  the  chemistry  model  were  also 
evaluated.  Comparison  of  solutions  assuming  fro¬ 
zen  and  kinetic  chemistry  approximations  indicates 
that  the  chemistry  model  influences  somewhat  the 
location  of  the  barrel  shock  reflection  point.  The 
other  features  of  the  plume  exhaust  flow  are  not 
drastically  affected  by  the  chemistry  assumption. 
However,  these  simulations  were  accomplished  at 
high-altitude  conditions,  where  the  entrainment  of 
the  atmospheric  gases  into  the  plume  shear  layer 
does  not  result  in  significant  afterburning.  It  is 
expected  that  chemistry  will  have  a  significant  influ¬ 
ence  on  solutions  at  lower  altitudes  where  plume/ 
atmospheric  afterburning  will  dominate  the  flow 


afterburning,  larger  effects  would  be  noted. 

It  should  also  be  noted  that  all  solutions 
assumed  constant  oxidizer/fuel  (0/F)  ratio  across 
the  nozzle  exit  plane.  When  non-ideal  engine 
effects  such  as  fuel-film  cooling  and  injector  imper¬ 
fections  are  accounted  for,  the  0/F  ratio  and,  con¬ 
sequently,  the  chemical  composition  of  the  exhaust 
products  can  vary  widely  across  the  nozzle  exit 
plane,  resulting  in  nonaxisymmetric  chemical  distri¬ 
butions.  Prediction  of  the  O/F  distribution  is  beyond 
the  scope  of  this  study  and  outside  of  the  current 
capabilities  of  the  GIFS  model. 

Three-Dimensional  Effect 

In  an  earlier  study  by  one  of  the  authors^  and 
others,®'^®  it  was  shown  that  three-dimensional 
effects  are  important  and  should  not  be  ignored  or 
oversimplified  in  modeling  efforts.  Previous  work 
focused  on  three-dimensional  plume  exhaust  phe¬ 
nomena  but  did  not  consider  the  upstream  influ¬ 
ence  of  the  missile  body,  missile  base  region,  or  the 
gas  generator  flow.  This  study  further  investigated 
the  three-dimensional  aspects  of  CFD  simulations, 
including  upstream  influences  generated  by  the 
Titan  II  missile  body,  base  region,  and  fuel-rich  gas 
generator  flow.  Static  pressure  and  Mach  number 
contours  resulting  from  these  three-dimensional 
calculations  initiated  at  the  missile  nose  are  shown 
in  Figs.  3  and  4,  respectively.  Figure  4  includes 
Mach  number  cross  sections  taken  at  the  100-m 


field.  The  chemical  kinetic  effect  strongly  influ¬ 
ences  the  calculation  in  the  immediate  missile 
base  region,  especially  if  the  fuel-rich  gas 
generator  flow  is  considered  in  the  simulation. 

The  influence  of  two-dimensional  versus 
one-dimensional  startline  conditions  was  also 
evaluated.  A  comparison  of  the  spatial  plume 
size  indicates  that  the  uniform  starting  line 
results  in  a  slightly  larger  plume  expansion 
region  due  to  the  increased  pressure  gradient 
between  the  nozzle  exit  and  the  free  stream. 


As  expected,  differences  between  the  two 
start  line  approximations  become  less  signifi¬ 
cant  as  the  flow  progresses  axially  down¬ 
stream.  Although  these  two  solutions  do  not 
show  a  significant  effect  of  exit  profile  shape 
on  the  plume  temperature,  it  is  expected  that 
in  other  cases,  especially  those  with  strong 


Fig.  3.  Pressure  contours  (3-D  calculations). 
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and  205-m  axial  locations  as  referenced  to  the  noz-  static  temperature  profiles  from  these  calculations 
zie  exit.  This  simulation  also  assumesf  finite,  rate-  are  compared  in  Fig.  6.  These  results  indicate  that 
controlled  chemistry  throughout  the  solution  the  plume  expansion  angles  are  larger  and  the 


domain,  including  the  fuel-rich  gas 
generator  flow  region.  The  plume 
expansion  shock  (barrel  shock)  for¬ 
mation  in  the  plume  near-field  and 
the  shock  reflection  visible  at  116 
meters  downstream  of  the  nozzle 
exit  plane  location  are  evident  in 
Figs.  3  and  4.  The  static  tempera¬ 
ture  Is  increased  Immediately 
downstream  of  the  reflection  point 
to  approximately  half  the  value  of 
the  total  temperature  of  the  flow. 
The  blunt  body  shock,  plume 
expansion  shock,  barrel  shock 
reflection  and  resulting  separation 
region  are  clearly  evident  In  these 
results.  Figures  5a  and  5b  are 
close-up  views  of  static  tempera¬ 
ture  contours  and  Mach  number 
vectors,  respectively,  in  the  missile 
base  flow  region.  The  heating  of 


*irrAH  H  m 


Fig.  4.  Mach  number  contours  for  3-D  reacting  flow  with  gas  generator. 


the  missile  base  surface  and  the  effect 
of  the  gas  generator  effluent  impinging 
on  the  nozzle  surfaces  are  clearly  evi¬ 
dent.  Further,  the  gas  generator  flow 
initially  expands  as  it  exits  the  nozzle, 
adjusting  to  the  ambient  pressure  con¬ 
dition,  but  further  downstream  in  the 
region  between  the  two  nozzles,  the 
gas  generator  flow  is  compressed 
because  of  the  area  change  created 
by  the  nozzle  expansion  skirts.  In  Fig. 
5b,  showing  the  Mach  vectors,  the 
interaction  of  the  gas  generator  flow 
with  the  base  region  and  the  nozzles  is 
apparent.  This  interaction  creates 
recirculation  of  the  hot  gases  into  the 
missile  base  region,  which  is  evident 
in  Fig.  5b. 

The  influence  of  the  missile  body 
and  base  regions  on  the  plume  flow 
field  is  assessed  by  comparing  solu¬ 
tions  obtained  with  and  without  the 


missile  body  and  base  region  a.  Static  temperature  contours 

included  in  the  simulation.  Static  tern-  pjg  5  jj^an  II  SLV  at  47.6  km  reacting  flow  calculation  near-field 

perature  contours  and  centerline  results. 


7 

American  Institute  of  Aeronautics  and  Astronautics 


overall  plume  width  is  larger  if  the  missile  body 
and  base  region  are  included  in  the  computational 
domain.  The  shock  reflection  point  was  located  12 
meters  farther  downstream  in  the  solution  includ¬ 
ing  the  body  and  base.  It  appears  that 
the  missile  body  effects  are  significant 
and  influence  the  plume  impingement 
shock  location,  the  inviscid  shock 
structure,  and  the  plume  shear  layer 
development. 

To  assess  the  influence  of  the  gas 
generator  effluents,  the  nose-to-tail 
solution  was  repeated  with  and  with¬ 
out  consideration  of  the  gas  generator 
flow  assuming  frozen  chemistry  condi¬ 
tions.  Comparison  of  these  results 
indicates  that  the  gas  generator  flow 
has  a  significant  influence  on  the  very 
near-field  portion  of  the  plume  flow 
field,  extending  to  approximately  1 
meter  downstream  of  the  nozzle  exit. 

The  effect  of  the  gas  generator  is  not 
noticeable  at  further  downstream  loca¬ 
tions.  To  assess  the  effect  of  the  gas 
generator.  Fig.  7  contrasts  Mach  num¬ 
ber  contours  and  radial  Mach  number 
profiles  taken  at  2-  and  3-m  axial  loca¬ 
tions  downstream  of  the  nozzle  exit 


plane  location  with  and  without  the  gas  generator 
flow.  The  gas  generator  flow  region  is  apparent 
near  the  centerline.  This  solution  assumed  a  frozen 
chemistry  condition  throughout  the  computationai 


Titan  tt  SLV  at  47.6  Km:  Reacting  Flow  Calculation 


A-.- 


syiii  v/ec  l.rgb 


b.  Mach  vectors 
Fig.  5.  Concluded. 


Fig.  6.  Static  temperature  contours  with  and  without  missiie  body. 
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i-in  7  Mflrh  niimhor  rnmnsiricnn  with  nflc  nonorfltnr  \/prciies  without  OaS  Generator. 

region.  Comparison  of  this  result  with  the  reacting  A  comparison  of  the  two-dimensional  axisym- 
gas  generator  results  indicates  that  in  order  to  metrio  solution  with  the  three-dimensiona!  twin- 
accurately  simulate  missile  base  heating  kinetic  nozzle  solution  was  accomplished  to  determine  the 
chemistry  is  paramount.  For  reacting  flow  condi-  impact  of  the  single  eguivalent  nozzle  assumption, 
tions,  it  is  also  expected  that  the  influence  of  the  This  simplification  is  often  used  to  approximate  and 
fuel-rich  gas  generator  flow  would  generate  strong  simplify  three-dimensional  geometries.  Figures  8a 
localized  combustion  zones  extending  to  further  and  8b  are  static  temperature  contours  and  provide 
downstream  positions.  a  comparison  of  the  axisymmstric  and  the  3-D 


a.  3-D  (jiuiiie  (top)  versus  single  equivalent  nozzle 


b.  3-D  (top)  versus  asisyniinetric  (bottom) 


(bottom) 


j-in  w  St3t!c  tsmpsrsturo  contour  compsrisons. 
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solutions.  In  both  figures,  the  3D  solution  is  shown  laminar  and  turbulence  stress  models  indicates 


33  tnS  uppSr  hsif  Of  ihS  COntOUr  plot  3nc!  tnS  3Xi“ 
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symmetric  solution  as  the  lower  portion  of  the  fig-  influenced  by  the  turbulence  model.  Laminar 


UI0.  I  11^1  OiO  oi^iiiiii^aiii  W4IIIC/I  i  me; 


lafSIl  ft  mlti 

oppi  v/>Mi  I  ictiivi  lo  Will  v^e^^ioue;  iuime;i  oi  iwve;t  ciiii~ 

two  solutions,  including  the  overall  plume  size,  tudes  where  turbulent  mixing  becomes  more  domi- 

asymmstric  flow  distributions,  and  the  location  of  nant.  In  general,  the  turbulent  approximation 

the  shock  reflection  points.  This  illustrates  that  a  resulted  In  slightly  increased  temperatures 

a  ttl  V%A^**Ia  AA  AI  ■  AA  AVAl  (I^AA  ^A  •  ■  A  A  •  A  flAtai  .*<J  A  W«  A  t  A 
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some  generalized  insight  concerning  overall,  nomi- 
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nal  cjuaiiiaiive  assessments;  I iuw«vtfi ,  lui  ueieneu  '-»■  ^/nesnn^uy  «oounipuwii  ve^i^uo 

studies  requiring  accurate  resolution  of  the  spatial  reacting)  did  not  influence  the  overall  plume  struc- 

character  of  mUitinOZZle  piUmc  nuwo,  a  unee-  luio  m  mi^  v^c*ow.  i  iwwnovoi ,  c*i  iv/wvyi-wimwvje; 

dimensional  calculation  is  required.  The  3-D  influ-  tions  where  afterburning  occurs,  the  chemistry 

r  .1  I  •  _  Xl.  ^  I  ..  .1.  .  lU _ —  —  aMaa+  VWA«r  1<«A  a!aAi4iAAA+  I  a  AAl^t‘4>AA  4|r^A  Al«Aa*AIA 

ence  is  further  rnagnilieo  wnen  me  uuuy/uclbe  cue  eiieui  may  we  oiymnwaiii.  Ill  awuluwii,  me  wiieimo- 
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reacting)  did  not  influence  the  overall  plume  struc- 


included  In  the  simulation. 


Conclusions 


A  modified,  three-dimensional  GIFS  computer 

AVAAk'AVA  BatAA  AAaI<A^  A%^4  A  A  A  S*  r  A  ll  *  4  A  ••  ^I^IA  A^ll^jBr 
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eiieei  iiiay  we  oi^iiinwcmi.  iii  owuiiiwii,  me  wiieiiiio- 

try  assumption  will  dramatically  influence  the  heat- 

^HAAA^AV  Ib'A  aI!aA  +  I  A  AA  IA  +1^  A  ArtlAAllA  l^AAA  ••AAIAA  AA^ 
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in  the  recirculation  and  impingement  regions 

beiween  rnuilipi^  no^zias. 


4.  Tha  jnucmOz.z.la  spacing  distanca  hcis  <a  sig- 


which  focused  on  evaluating  the  three-dimensional  nificant  impact  on  the  barrel  shock  formation  and 


effects  of  the  missile  body*,  miSSlIe  wSSe  TeQICn^  ond 


raf lection  location,  pluma/pluina  impinganiant 


the  gas  generator  on  the  resulting  plume  exhaust  shock  location,  and  shear  layer  development. 
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Instances  where  these  influences  may  impact 


resulting  from  a  Titan  11  propulsion  system  operat-  flow-field  predictions  occur  when  nozzle  gim- 
ing  at  high-altitude  (47.6  km)  flight  conditions  was  balling  cycles  are  considered, 
parametrically  simulated  in  this  numerical  invest!- 

aa4SaA  D*>AB|!ai>A  a4|  I^IAA  •  A«  « AA  +  iA  a4a^  A  aM  A  A+A  C  A  A  A  A>«  A  A  A  A  A  a4  ^t^A  4»«(A  ^IAAAAaSaAaI  AaW 

gatlCn.  I  r0VtOuS  SiUGI0S  mveSu^aiew  me  eiiewie  w.  r\  wwi  i  i^cii  lowi  i  wi  me  iww~wii  i  lei  leiwMoi  oiiu 

resulting  from  startline  nonuniformity,  intranozzle  three-dimensional  cases  indicates  that  the  three- 
spacing,  turbulence,  and  finite-rate  chemistry  for  dimensional  effects  are  important  in  the  near-field 
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the  same  Titan  II  twin -nozzle/plume  propulsion  plume  and  diminish  as  the  axial  distance  extends 
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this  study.  Analysis  of  all  the  results,  including  the  The  single,  equivalent  nozzle  approach  should  not 
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teristics  where  three-dimensional  effects  are  pre¬ 


file  versus  one-dimensional  uniform,  had  a  signifi-  tures  and  spatial  information  are  required  as  part 
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culation.  The  difference  between  the  two  startline 
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distance  downstream  of  the  nozzle  exit  plane  was  the  barrel  shock  reflection  location,  plume/piume 
increased.  This  significance  of  the  observed  sensi-  impingement  shook  location,  and  development  of 
tivity  needs  additional  investigation  to  elucidate  the  the  shear  layer  region.  The  shock  layer  initiated 
effect  of  the  starting  ccnditions  on  the  plume  flow  near  the  missile  ncse  increased  the  temperature 

field,  especially  in  the  presence  of  strong  afterburn-  surrounding  the  missile  body.  This  temperature 
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the  shear  layer  region.  The  shock  layer  initiated 
noGr  tnG  rrii33i!0  nc30  irior0330cl  thi0  t0rTip0r3tur0 
surrounding  the  missile  body.  This  temperature 
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which  are  expected  to  dominate  lower  altitude  sim-  field  domain. 
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7.  A  comparison  of  solutions  contrasting  gas 
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flow  has  a  pronounced  influence  in  the  immediate 
plume  near-field  region  and  diminishes  beyond 
approximately  1  meter  downstream  of  the  nozzle 
exit  location.  The  extent  of  the  effect  is  expected 
to  be  based  on  the  gas  generator  mass  fiow.  The 
gas  generator  flow  is  expected  to  contribute  to  the 
missile  base  heating  and  nozzle  wall  impingement 
heating,  resulting  in  recirculation  of  hot  gas  in  the 
region  between  the  two  nozzles.  However,  the  cur¬ 
rent  frozen  chemistry  solution  did  not  adequately 
simulate  the  heating  effect. 

8.  These  analyses  indicate  that  two-dimen¬ 
sional,  axisymmetric,  nonequilibrium  analysis  tools 
can  provide  generai  insight  concerning  general¬ 
ized  qualitative  assessments  of  multiple  plume 
flow-field  phenomena.  However,  for  detailed  stud¬ 
ies  of  complex  flow-field  phenomena  and  spatial 
generalizations,  a  more  sophisticated  three-dimen¬ 
sional  calculation  is  required. 

Furthermore,  this  study  is  intended  to  assess, 
understand,  and  quantify  propulsion  and  plume 
physical  phenomena  in  order  to  identify  where 
simplified  models  can  be  used  without  introducing 
significant  errors  in  conclusions  that  might  be 
deduced  from  analysis  of  the  flow-field  simulation. 
Toward  this  objective,  this  study  explored  three- 
dimensional  vehicle/base  interactions  with  the 
plume  flow.  Also,  to  promote  optimal  utilization  of 
computer  resources,  hybrid  Navier-Stokes/ 
Parabolized  Navier-Stokes  methodologies  have 
proven  to  be  very  efficient. 
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